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New  Results  on  Modal  Participation  Factors: 
Revealing  a  Previously  Unknown  Dichotomy 

Wael  A.  Hashlamoun,  Munther  A.  Hassouneh  and  Eyad  H.  Abed 


Abstract — This  paper  presents  a  new  fundamental  approach 
to  modal  participation  analysis  of  linear  time-invariant  systems, 
leading  to  new  insights  and  new  formulas  for  modal  participation 
factors.  Modal  participation  factors  were  introduced  over  a  quar¬ 
ter  century  ago  as  a  way  of  measuring  the  relative  participation  of 
modes  in  states,  and  of  states  in  modes,  for  linear  time-invariant 
systems.  Participation  factors  have  proved  their  usefulness  in  the 
field  of  electric  power  systems  and  in  other  applications.  However, 
in  the  current  understanding,  it  is  routinely  taken  for  granted 
that  the  measure  of  participation  of  modes  in  states  is  identical 
to  that  for  participation  of  states  in  modes.  Here,  a  new  analysis 
using  averaging  over  an  uncertain  set  of  system  initial  conditions 
yields  the  conclusion  that  these  quantities  (participation  of  modes 
in  states  and  participation  of  states  in  modes)  should  not  be 
viewed  as  interchangeable.  In  fact,  it  is  proposed  that  a  new 
definition  and  calculation  replace  the  existing  ones  for  state  in  mode 
participation  factors,  while  the  previously  existing  participation 
factors  definition  and  formula  should  be  retained  but  viewed 
only  in  the  sense  of  mode  in  state  participation  factors.  Several 
examples  are  used  to  illustrate  the  issues  addressed  and  the 
results  obtained. 

Index  Terms — Participation  factors,  modal  participation  fac¬ 
tors,  modal  analysis,  linear  systems,  stability,  control  systems. 


I.  Introduction 

This  paper  presents  new  concepts,  results,  and  formulas 
in  the  subject  of  modal  participation  analysis  of  linear  time- 
invariant  systems.  This  topic  is  an  important  component  of 
the  Selective  Modal  Analysis  (SMA)  framework  introduced 
by  Perez- Arriaga,  Verghese  and  Schweppe  [7],  [12]  in  the 
early  1980s.  A  main  construct  in  SMA  is  the  concept  of 
modal  participation  factors  (or  simply  participation  factors). 
Participation  factors  are  scalars  intended  to  measure  the  rel¬ 
ative  contribution  of  system  modes  to  system  states,  and  of 
system  states  to  system  modes,  for  linear  systems.  The  work  of 
these  authors  has  had  a  major  impact  especially  in  applications 
to  electric  power  systems,  where  participation  factors  as  they 
were  originally  introduced  have  become  a  routine  tool  for  the 
practitioner  and  researcher  alike. 

Since  their  introduction,  participation  factors  have  been 
employed  widely  in  electric  power  systems  and  other  ap¬ 
plications.  They  have  been  used  for  stability  analysis,  order 
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reduction,  sensor  and  actuator  placement,  and  coherency  and 
clustering  studies  (e.g.,  [7],  [12],  [8],  [2],  [5],  [3],  [9]).  Several 
researchers  have  also  considered  alternate  ways  of  viewing 
modal  participation  factors  (e.g.,  [11],  [4],  [10]). 

We  study  linear  time-invariant  continuous-time  systems 

x  =  Ax{t)  (1) 

where  x  £  R"  and  A  is  a  real  n  x  n  matrix.  We  make  the 
blanket  assumption  that  A  has  a  set  of  n  distinct  eigenvalues 
The  solution  of  (1)  then  takes  the  form  of  a 
sum  of  modal  components: 

x{t)  =  £ex,V  (2) 

i=i 

where  the  c'  are  constant  vectors  determined  by  the  initial 
condition  x°  and  by  the  right  and  left  eigenvectors  of  A. 

In  their  study  of  modal  participation  for  the  system  (1),  the 
authors  of  [7],  [12]  selected  particular  initial  conditions  and 
introduced  definitions  motivated  by  the  calculation  of  relative 
state  and  mode  contributions  using  those  initial  conditions. 
In  this  paper,  we  take  a  different  approach,  building  on  our 
previous  work  [1],  in  which  definitions  of  modal  participation 
factors  are  formulated  by  averaging  relative  contributions  of 
modes  in  states  and  states  in  modes  over  an  uncertain  set  of 
initial  conditions.  In  this  approach,  we  consider  initial  condi¬ 
tions  to  be  unknown,  and  we  take  the  view  that  performing 
some  sort  of  average  over  all  possible  initial  conditions  should 
give  a  more  reliable  result  than  focusing  attention  on  one 
particular  possible  initial  condition.  The  uncertainty  in  initial 
condition  can  be  taken  as  set-theoretic  (unknown  but  bounded) 
or  probabilistic.  We  took  the  same  basic  approach  in  our 
paper  [1],  but  later  found  a  subtle  error  in  the  calculation  in 
that  paper  for  the  case  of  state-in-mode  participation  factors. 
Upon  realizing  this  subtle  error,  we  embarked  on  the  present 
research,  in  which  we  find  a  previously  unnoticed  dichotomy 
in  the  two  basic  types  of  modal  participation  factors. 

The  main  contribution  of  this  paper  is  to  reveal  this  pre¬ 
viously  unknown  dichotomy  in  modal  participation  analysis. 
To  wit,  although  the  definitions  obtained  in  [7],  [12],  and 
which  have  been  in  wide  use  since  their  introduction,  give 
identical  values  for  measures  of  participation  of  modes  in 
states  and  for  participation  of  states  in  modes,  these  are 
in  fact  better  viewed  as  fundamentally  different,  and  should 
be  calculated  using  two  distinct  formulas.  Summarizing,  the 
main  contribution  of  this  paper  is  as  follows:  we  propose 
replacing  the  existing  definition  of  participation  factors  with 
two  separate  definitions  that  yield  distinct  numerical  values 
for  participation  of  modes  in  states  and  for  participation  of 
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states  in  modes.  In  this  paper,  the  currently  used  participa¬ 
tion  factors  measuring  participation  of  states  in  modes  are 
replaced  with  a  new  first-principles  definition,  a  particular 
instance  of  which  is  an  explicit  formula  given  in  Section 
V.  In  addition,  we  show  that  our  formula  for  participation 
factors  measuring  participation  of  modes  in  states  agrees 
with  the  commonly  used  participation  factors  formula  under 
reasonable  assumptions  on  the  allowed  uncertainty  in  the 
system  initial  conditions.  Thus,  a  dichotomy  is  proposed  in 
the  calculation  of  participation  factors. 

The  paper  proceeds  as  follows.  In  Section  II,  the  original 
definitions  of  modal  participation  factors  are  recalled  from  [7], 
[12],  In  Section  III,  basic  examples  are  used  to  illustrate 
the  need  for  an  approach  that  yields  distinct  formulas  for 
measuring  the  two  main  types  of  modal  participation:  par¬ 
ticipation  of  modes  in  states  and  participation  of  states  in 
modes.  In  Section  IV,  the  approach  we  introduced  to  this 
topic  in  [1]  is  recalled  and  discussed  in  light  of  the  objectives 
of  the  present  paper.  The  discussion  makes  clear  that  this 
approach,  based  on  defining  modal  participation  measures  by 
averaging  over  an  uncertain  set  of  initial  conditions,  readily 
yields  the  original  definition  for  participation  of  modes  in 
states,  but  does  not  easily  yield  a  simple  closed-form  expres¬ 
sion  for  measuring  participation  of  states  in  modes.  Next,  in 
Section  V,  a  candidate  closed-form  formula  is  obtained  for 
modal  participation  factors  that  measure  the  participation  of 
states  in  modes;  this  is  achieved  by  careful  evaluation  of  the 
general  averaging  formula  from  Section  IV  under  a  simplifying 
assumption  on  the  initial  condition  uncertainty.  It  is  important 
to  note  that  to  obtain  this  simple  formula,  a  specific  form  is 
assumed  for  the  uncertainty  in  the  system  initial  condition; 
other  assumptions  on  the  initial  condition  uncertainty  would 
not  lead  to  the  same  formula  or  any  readily  useable  expression. 
The  derived  formula  is  proposed  since  it  reflects  the  effect  of 
initial  condition  uncertainty  and  can  be  derived  analytically.  In 
Section  VI,  a  mechanical  system  example  is  used  to  illustrate 
the  usefulness  of  the  explicit  formula  for  state  in  mode  par¬ 
ticipation  factors.  In  Section  VII,  an  additional  result  is  given 
that  relates  only  to  mode  in  state  participation  factors.  This 
result  expands  the  initial  condition  uncertainty  assumptions 
under  which  the  traditional  participation  factors  formula  can 
be  shown  to  accurately  measure  mode  in  state  participation 
using  the  averaging  formulation.  Concluding  remarks  and 
suggestions  for  future  work  are  collected  in  Section  VIII. 

II.  Original  Definitions  of  Modal  Participation 
Factors 

In  this  section,  the  original  definitions  of  modal  participa¬ 
tion  factors  are  recalled  from  [7],  [12],  Consider  the  linear 
system  (1),  repeated  here  for  convenience: 

x  =  Ax(t)  (3) 

where  x  £  R",  and  A  is  a  real  n  x  n  matrix.  The  authors 
of  [7],  [12]  also  make  the  blanket  assumption  that  A  has 
n  distinct  eigenvalues  Let  (rl  ,r2, . . .  ,rn)  be 

right  eigenvectors  of  the  matrix  A  associated  with  the  eigen¬ 
values  (A-i , A-2, . . - , A.„),  respectively.  Let  (I1,/2,...,/")  denote 


left  (row)  eigenvectors  of  the  matrix  A  associated  with  the 
eigenvalues  (A-i,A-2,  • . .  ,A.„),  respectively.  The  right  and  left 
eigenvectors  are  taken  to  satisfy  the  normalization  [6] 

/V  =  8y  (4) 

where  8 y  is  the  Kronecker  delta: 


The  solution  to  (3)  starting  from  an  initial  condition  x(0)  =  x° 
is 

x(t)  =  eA,x°  (5) 

Since  the  eigenvalues  of  A  are  distinct,  A  is  similar  to  a 
diagonal  matrix.  Using  this,  (5)  can  be  rewritten  in  the  form 

=  !(/'*  V‘V.  (6) 

(=i 

From  (6),  xt ft)  is  given  by 

xk(t)  =  i(lixy^ri.  (7) 

i=  1 

A.  Relative  participation  of  the  i-th  mode  in  the  k-th  state 

To  determine  the  relative  participation  of  the  i-th  mode 
in  the  k-th  state,  the  authors  of  [7],  [12]  select  an  initial 
condition  x°  =  ek,  the  unit  vector  along  the  k-th  coordinate 
axis.  As  seen  next,  this  choice  is  convenient  in  that  it  results 
in  a  simple  formula  for  mode-in-state  participation  factors.  We 
note  that  the  derived  formula  for  mode-in-state  participation 
factors  agrees  with  that  obtained  using  an  uncertain  initial 
condition  under  general  assumptions,  as  demonstrated  in  [1] 
and  in  Section  IV  below.  With  this  choice  of  x°,  the  evolution 
of  the  k-th  state  becomes 

**(o  =  £(/&)«* 

i=i 

=:  (8) 

1=1 

The  quantities 

Pki  :=  l‘A  (9) 

are  found  to  be  unit-independent,  and  are  taken  in  [7],  [12]  as 
measures  of  the  relative  participation  of  the  i-th  mode  in  the 
k-th  state;  p *,■  is  defined  in  [7],  [12]  as  the  participation  factor 
for  the  i-th  mode  in  the  k-th  state. 

B.  Relative  participation  of  the  k-th  state  in  the  i-th  mode 
The  relative  participation  of  the  k-th  state  in  the  i-th  mode 

is  studied  in  [7],  [12]  by  first  applying  the  similarity  transfor¬ 
mation 

z  :=  V_1x  (10) 

to  system  (3),  where  V  is  the  matrix  of  right  eigenvectors  of 
A: 


V  =  [r1  r2  ■■■  rn] 


(11) 
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and  V  1  is  the  matrix  of  left  eigenvectors  of  A: 

[  l 1  1 

v~l  =  .  . 

.  r  . 

Then  z  obeys  the  dynamics 

z(t)  =  V~lAVz(t) 

=  A z(t), 


(12) 


(13) 


where  A  :=  diag(Ai,  X2,  ■  ■  ■ ,  A„),  with  initial  condition  z°  := 
y-‘x°.  This  implies  that  the  evolution  of  the  new  state  vector 
components  zu  t  =  1 .....  n  is  given  by 


Zi(t) 


zy-> 

lix°ex>t 


n 


I  (/**?) 


Ai? 


(14) 


For  a  real  eigenvalue  A;,  clearly  Zi(t)  represents  the  evolution 
of  the  associated  mode.  If  A;  is  not  real,  then  the  associated 
mode  is  sometimes  taken  to  be  Zi(t),  but  can  also  be  taken  as 
the  combination  of  Zi(t)  and  its  complex  conjugate  z* ( t ),  which 
reflects  the  influence  of  the  eigenvalue  X*.  In  the  latter  ap¬ 
proach,  we  view  A,  and  A*  as  representing  the  same  “complex 
frequency.”  In  the  past,  the  former  convention  was  used  in  most 
publications.  In  this  paper,  we  allow  both  interpretations,  but 
we  will  find  it  convenient  to  use  the  latter  point  of  view  when 
deriving  a  new  state-in-mode  participation  factors  formula  for 
the  case  of  complex  eigenvalues. 

In  order  to  determine  the  relative  participation  of  the  k- 
th  state  in  the  ;-th  mode,  the  authors  of  [7],  [12]  select  an 
initial  condition  x°  =  r\  the  right  eigenvector  associated  with 
A,-.  As  seen  next,  this  choice  is  convenient  in  that  it  results  in 
a  simple  formula  for  state-in-mode  participation  factors.  We 
will  revisit  this  later  using  an  uncertain  initial  condition,  and 
obtain  a  different  result.  With  this  choice  of  initial  condition, 
the  evolution  of  the  /-th  mode  becomes 


Zi(t) 


/VeV 


k'k 


k=  1 


Efti 

k=\ 


At 


At 


(15) 


Based  on  (15),  the  authors  of  [7],  [12]  propose  the  formula 


Pki  =  I'A 


(16) 


as  a  measure  of  the  relative  participation  of  the  A-th  state  in 
the  ;-th  mode. 

Note  that  (9),  (16)  provide  identical  formulas  for  participa¬ 
tion  of  modes  in  states  and  participation  of  states  in  modes, 
respectively.  For  this  reason,  the  same  notation  pia  was  used 
for  both  types  of  participation  factors  until  now. 


III.  Motivating  Examples  Showing  Inadequacy  of 
Participation  Factors  Formula  as  a  Measure  of 
State  in  Mode  Participation 

In  this  section,  by  way  of  motivation  for  the  subsequent 
analysis,  two  examples  are  given  that  show  the  need  for  a  new 
definition  and  a  new  formula  for  state  in  mode  participation 
factors. 


Example  1  Consider  the  two-dimensional  system 


X\ 

a 

b 

XI 

.  *2 

0 

d 

.  X2 

A 


where  a ,  b  and  d  are  constants  with  a  /  d.  The  eigenvalues 
of  A  are  Aj  =  a  and  X2  =  d.  The  right  eigenvectors  associated 
with  Aj  and  At  are 


r1  = 


and  rl  = 


1 

d—a 

b 


respectively.  The  left  eigenvectors  associated  with  Aj  and  At 
and  satisfying  the  normalization  (4)  are 


/' 


1 


b 

a—d 


and  l 2 


0 


-b 

a—d 


respectively. 

Before  calculating  the  participation  factors  measuring  the 
influence  of  states  x\  and  at  in  mode  1 }  we  write  the  evolution 
of  mode  1  explicitly.  Using  (14),  we  have 

zi(f)  =  llx°eXlt 


a  —  d 

=  (4  +  ^?)^.  (17) 

Note  that  the  evolution  of  mode  1  is  influenced  by  both  Xj 
and  Xj,  with  the  relative  degree  of  influence  depending  on  the 
values  of  the  system  parameters  a,b  and  d. 

Calculating  the  participation  factors  using  the  original  def¬ 
inition  as  recalled  in  the  foregoing  section,  we  find  the 
participation  factor  for  state  x\  in  mode  1  is  p\\  =  l\r\  =  1, 
while  the  participation  factor  for  state  X2  in  mode  1  is 
P2\  =  l\r\  =  0.  Thus,  the  original  definition  of  participation 
factors  for  state  in  mode  participation  indicates  that  state  at 
has  much  smaller  (even  zero)  influence  on  mode  1  compared 
to  the  influence  coming  from  state  xi ,  regardless  of  the  values 
of  system  parameters  a,b  and  d.  This  is  in  stark  contradiction 
to  what  we  observed  using  the  explicit  formula  (17),  and  begs 
for  a  re-examination  of  the  basic  formula  for  state-in-mode 
participation  factors.  ■ 


M1 


^  [  or  simplicity,  we  use  the  terminology  'mode  f  in  place  of  'the  mode 
associated  with  eigenvalue  Xj.’ 
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Example  2  Consider  the  two-dimensional  system 


x\ 

1 

1 

XI 

.  *2 

-d 

-d 

.  X2  . 

A 


where  d  ^  1  is  a  constant.  The  eigenvalues  of  A  are  Aj  =  0 
and  A-2  =  1  —d.  The  right  eigenvectors  associated  with  Ai  and 
A.  2  are 


r 


l 


1 

-1 


and  r2 


1 

-d 


respectively.  The  left  eigenvectors  associated  with  A]  and  A2 
and  satisfying  the  normalization  (4)  are 


/' 


r  _=d_  _^t_ 

.  1  -d  1  -d 


and  / 2 


1 

1  -d 


1 

1  —d 


respectively.  Denote  by  V  the  matrix  of  right  eigenvectors  of 
A: 


V  =  [r'  r2} 


1  1 

-1  -d 


From  the  normalization  condition  (4),  we  can  immediately 
write 


r  il  ] 

r  -a  -1  1 

1  _ 

1 -d  1  -d 

12  _ 

1  1 

.  1  -d  1  -d  . 

The  evolution  of  the  modes  can  be  obtained  using  the 
diagonalizing  transformation  z  ■=  V  lx  as  was  done  in  (10)- 
(13).  The  system  modes  are  found  to  be 


zi(t) 

■  [lx0 eXlt  - 

.  Z2(0  . 

l2x°ex2t 

(i  .. 
T^C 


-d  r0 _ 1 

dX  1 
.0 


■0)  P~X it 


1 


1  -<#-*2) e 


*2)  e 


(18) 


Based  on  the  original  definition  of  participation  factors,  the 
participation  factor  for  state  x\  in  mode  2  is  p\2  =  r\l 2  = 
and  the  participation  factor  for  state  X2  in  mode  2  is  P22  = 
r2^2  =  T^Fr  Clearly,  in  general  p  12  ^  P22-  However,  from  (18) 
we  have  the  equation 

Z2{t)  =  ^{x\+x§ex*  (19) 

for  the  second  mode  22(f),  from  which  we  observe  that  state 
x\  and  state  X2  participate  equally  in  mode  2  since  22(f) 
depends  on  the  initial  condition  x°  through  the  sum  x®  +x®. 
Again,  we  find  that  the  state-in-mode  participation  factors  as 
commonly  calculated  yield  conclusions  that  are  very  much 
at  odds  with  what  one  might  consider  reasonable  based  on 
explicit  calculation  of  the  evolution  of  system  modes  as  they 
depend  on  initial  conditions  of  the  state  variables.  ■ 


IV.  Initial  Condition  Uncertainty  Approach: 
Application  to  Derivation  of  Mode-in-State 
Participation  Factors 


For  systems  operating  near  equilibrium,  it  is  often  reason¬ 
able  to  view  the  system  initial  condition  as  being  an  uncertain 
vector  in  the  vicinity  of  the  system  equilibrium  point.  In  this 
paper,  and  in  the  authors’  previous  work  [1],  we  approach 
the  problem  of  measuring  modal  participation  by  averaging 
relative  contributions  over  an  uncertain  set  of  initial  conditions. 
In  this  section,  we  summarize  this  approach  as  it  applies  to 
the  definition  and  calculation  of  mode-in-state  participation 
factors.  We  carry  the  approach  through  to  its  conclusion 
for  this  problem,  obtaining  an  explicit  formula  for  mode- 
in-state  participation  factors.  As  mentioned  above,  the  final 
formula  we  obtain  using  this  approach  agrees  in  this  case 
with  the  previously  existing  expression  p^.  However,  in  the 
next  section,  such  a  happy  coincidence  will  not  occur  for  the 
more  delicate  situation  of  defining  and  calculating  state-in- 
mode  participation  factors. 

Next,  we  recall  from  our  previous  work  [1]  a  basic  definition 
of  relative  participation  of  a  mode  in  a  state.  This  definition 
involves  taking  an  average  over  system  initial  conditions  of  a 
measure  of  the  relative  influence  of  a  particular  system  mode 
on  a  system  state.  The  initial  condition  uncertainty  can  be 
taken  as  set-theoretic  or  probabilistic.  In  the  set-theoretic  for¬ 
mulation,  the  participation  factor  measuring  relative  influence 
of  the  mode  associated  with  A,-  on  state  Xk  can  be  defined  as 


Pki 


avg 
x°  e  § 


(l'x°)r'k 


(20) 


whenever  this  quantity  exists  (however,  see  Remark  1  below 
for  another  possible  definition  for  the  case  of  complex  A,). 
Here,  x®  =  L/Li  (l'x°)r'k  is  the  value  of  x^t)  at  t  =  0,  and 
“avgvoG§”  is  an  operator  that  computes  the  average  of  a 
function  over  a  set  §  C  R"  (representing  the  set  of  possible 
values  of  the  initial  condition  x°).  We  assume  that  the  initial 
condition  uncertainty  set  §  is  symmetric  with  respect  to  each 
of  the  hyperplanes  {v>  =  0},  k  =  1 

In  the  definition  in  [1]  that  starts  with  a  probabilistic 
description  of  the  uncertainty  in  the  initial  condition  x°,  the 
average  in  (20)  is  replaced  by  a  mathematical  expectation. 
The  general  formula  for  the  participation  factor  pu  measuring 
participation  of  mode  i  in  state  Xk  becomes 


Pki 


E 


J  (i‘x°K  \ 

I  4  J 


(21) 


The  inadequacy  of  the  original  state-in-mode  participation 
factors  formula  has  been  demonstrated  in  the  two  examples 
above.  This  motivates  the  need  for  a  new  formula  that  better 
assesses  the  influence  of  system  states  on  system  modes. 


where  the  expectation  is  evaluated  using  some  assumed  joint 
probability  density  function  f(x°)  for  the  initial  condition 
uncertainty  (of  course,  this  definition  applies  only  when  the 
expectation  exists). 
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Expanding  the  inner  product  term  in  (21),  we  find 


Pki 


E 

E 


7=1  Xk  I 


IMH 
l  4 


+  E 


iiA+  t  l‘AE 

7=i,  m 


E 

7=1,  jAk 


AfM  \ 

4  ) 


(22) 


The  second  term  in  (22)  vanishes  when  the  components  of 
the  initial  condition  vector  x^x®,  . . .  ,xjj  are  independent  with 
zero  mean  [1],  Therefore,  under  the  assumption  that  the  initial 
condition  components  Xj  ,x?, . . . ,  xjj  are  independent  with  zero 
mean,  the  participation  of  the  ;-th  mode  in  the  A-th  state  is 
given  by  the  same  expression  originally  introduced  by  Perez - 
Arriaga,  Verghese  and  Schweppe  [7],  [12]: 


ii  i 

Pki  =  lkrk- 


(23) 


This  result  can  also  be  obtained  using  the  set-theoretic  aver¬ 
aging  formula  (20)  [1], 

Remark  1:  ( Alternate  Definition  of  Mode-in-State  Partici¬ 
pation  Factor  for  a  Complex  Mode )  For  a  complex  eigen¬ 
value  A the  associated  “mode”  is  taken  above  as  the  term 
containing  ek'1  in  the  system  response  (2).  However,  we  can 
alternately  view  this  mode  as  consisting  of  the  combined 
contributions  from  A,  and  its  complex  conjugate  eigenvalue  A*. 
This  viewpoint  is  easily  seen  to  lead,  under  the  same  symmetry 
hypotheses  as  above,  to  the  following  alternate  expression  for 
the  participation  factor  of  the  mode  associated  with  A;  and  A* 
in  state  xk: 

Pki  =  2Re{Zi4}.  (24) 


V.  New  Definition  of  Participation  Factors 
Measuring  Participation  of  States  in  Modes 

In  this  section,  a  new  definition  and  calculation  are  given 
for  participation  factors  measuring  contribution  of  states  in 
modes.  The  probabilistic  approach  presented  in  the  previous 
section  is  used,  where  the  initial  condition  is  assumed  to  satisfy 
a  joint  probability  density  function.  In  order  to  obtain  an 
explicit  formula  from  the  new  general  definition  of  state-in- 
mode  participation  factors,  we  find  that  it  is  necessary  to  make 
an  assumption  on  the  probability  distribution  of  the  initial  con¬ 
dition  which  is  more  constraining  than  what  was  needed  in  the 
analysis  above  for  mode-in-state  participation  factors.  Thus, 
the  explicit  formula  derived  in  this  section  should  be  viewed  in 
the  pragmatic  sense  that  it  provides  an  easy  to  use  expression 
that  reflects  initial  condition  uncertainty.  Other  assumed  forms 
of  uncertainty  may  not  lead  to  explicit  formulas,  although  a 
formula  requiring  numerical  evaluation  of  integrals  can  always 
be  obtained  from  the  definition.  The  explicit  formula  obtained 
here  differs  from  the  single  formula  (16)  that  is  currently 
used  to  measure  both  state-in-mode  participation  and  mode- 
in-state  participation,  while  the  currently  used  formula  (16) 
is  retained  here  as  a  measure  of  mode-in-state  participation 
(noting  that  the  alternate  formula  (24)  can  also  be  used  for 
the  case  of  a  complex  mode).  This  dichotomy  represents  a 


significant  departure  from  current  practice.  We  will  also  use 
the  new  formula  to  revisit  the  examples  of  Section  III. 

Consider  the  general  linear  time-invariant  continuous-time 
system  given  in  (3),  repeated  here  for  convenience: 

x  =  Ax(t)  (25) 


The  evolution  of  Zi{t),  i  =  1 was  obtained  in  Section  II 
as 


Zi{t)  =  eXi,l‘x° 

=  £$*?)■  (26) 
7=1 


This  equation  shows  the  contribution  of  each  component  x®, 
j  =  1, . . .  ,n  of  the  initial  state  x°  to  Zi(t).  Recall  also  that  for  the 
case  of  a  real  eigenvalue  A/,  Zi{t)  is  identically  the  i-th  mode, 
while,  for  a  complex  eigenvalue  A,,  the  associated  mode  can  be 
taken  as  Ziit)  or  as  the  combination  of  zft)  and  its  conjugate: 
Zi(t)+Z*(t)  =2Re  {z,(0}-  The  following  general  definition  of 
state-in-mode  participation  factors  is  obtained  by  averaging  the 
relative  contribution  of  x?  in  the  i-th  mode  and  evaluating  the 
result  at  t  =  0.  In  this  definition,  we  take  the  mode  associated 
with  a  complex  eigenvalue  as  2Re  {zft)},  i.e.,  the  combination 
of  modal  components  due  to  the  eigenvalue  and  its  conjugate. 
Had  we  decided  to  view  the  mode  associated  with  a  complex 
eigenvalue  A;  as  z,(t)  alone,  we  would  use  the  first  expression 
in  the  definition  below  for  both  the  case  of  a  real  and  a 
complex  eigenvalue.  However,  the  derivation  following  the 
basic  definition  below  of  a  simple  final  formula  would  become 
unwieldy  for  the  complex  eigenvalue  case. 

Definition  1:  For  a  linear  time-invariant  continuous-time 
system  (25),  the  participation  factor  for  the  A-th  state  in  the 
/-th  mode  is 


Ttkj 


E 

E 


JMH 


if  A,-  is  real 
if  A,  is  complex 


(27) 


whenever  the  expectation  exists. 

Note  that  in  (27),  the  notation  z®  means  zft  =  0)  =  l‘x°  and 
the  asterisk  denotes  complex  conjugation.  Also,  analogous  to 
the  approach  in  Section  IV  and  the  original  work  [7],  [12],  the 
quantities  being  evaluated  represent  the  contribution  of  state  xk 
to  a  mode  divided  by  the  total  mode  evaluated  at  time  7=0 
(however,  see  the  Conclusions  section  about  the  possibility 
of  measuring  modal  participation  effects  over  time).  Note 
also  that  Definition  1  always  yields  real- valued  participation 
factors. 

Unfortunately,  even  under  an  assumption  such  as  symmetry 
of  the  initial  condition  uncertainty,  there  is  no  single  closed- 
form  expression  for  the  state  in  mode  participation  factors 
Kkk  To  obtain  a  simple  closed-form  expression  for  the  state 
in  mode  participation  factors  nki  using  (27),  we  need  to 
find  an  assumption  on  the  probability  density  function  /(x°) 
governing  the  uncertainty  in  the  initial  condition  x°  that  allows 
us  to  explicitly  evaluate  the  integrals  inherent  in  the  definition. 

In  the  remainder  of  this  section,  we  assume  that  the  units 
of  the  state  variables  have  been  scaled  to  ensure  that  the 
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probability  density  function  /(x°)  is  such  that  the  components 
Xj,X2,...,x®  are  jointly  uniformly  distributed  over  the  unit 
sphere  in  Rn  centered  at  the  origin: 


f(x°) 


k  f|x°||<l 

0  otherwise 


(28) 


(This  is  the  same  as  assuming  a  uniform  distribution  in 
an  ellipsoid  that  is  centered  at  the  origin  and  symmetric 
with  respect  to  the  coordinate  hyperplanes  in  the  original 
state  variable  units,  a  physically  palatable  assumption  and 
independent  of  units  by  construction.)  The  constant  k  is  chosen 
to  ensure  the  normalization 

L<Jix°)dx°=1'  <29) 


The  value  of  the  constant  k  can  be  determined  by  evaluating 
the  integral  in  (29)  using  /(x°)  given  in  (28): 


f  f(x°)dx°=  f  kdx°ldx%...dx°n=kVn  =  l  (30) 

^ikon<i  4°h<i 


where  Vn  is  the  volume  of  the  unit  sphere  in  R" .  The  constant 
k  is  then  given  by 


k  = 


1 

Vn 


(31) 


The  following  Lemma  will  be  used  below. 

Lemma  1:  For  vectors  a  £  C",  b  £  Rn  with  b  ^  0  we  have 


T 

a  x 

Wid'x 


iTb 


bTb 


V„ 


(32) 


where  d„x  denotes  the  differential  volume  element 
dx\dx2  ■  ■  dxn,  and  Vn  is  the  volume  of  a  unit  sphere  in 
R"  which  is  given  by 


V„  = 


2, 

", 


n  =  1 
n  =  2 
n  >  3 


(33) 


Proof:  The  proof  below  is  for  the  case  a, b  £  Rn,  b  /  0.  The 
case  a  £  Cn  and  b  £  R"  follows  by  linearity.  Consider  the 
transformation  x  =  Qy  where  the  matrix  Q  is  chosen  to  be  an 
orthogonal  matrix  (i.e.,  Q =QT)  with  first  column  b1  :=  pjj 
(i.e.,  Q  =  [bl  Q\}).  The  remaining  columns  of  Q\  are  chosen 
to  be  orthogonal  to  bl,  i.e., 

C bl)TQi=0 .  (34) 

With  this  transformation,  since  Q  1  =  Q1,  we  have  that 
|  |x|  |2  =  xrx  =  (Qy)T Qy  =  yTQTQy  =  yTy  =  |  |y 1 12.  Also,  since 
Q  is  orthogonal,  det(Q)  =  1  and  d„x  =  dny. 

Writing  the  vector  y  as 


y  = 


y  i 
/ 

y 


with  yi  £R  and  y 
as 


l  Pd-X  =  / 

J \\x\\<\blx  J\\y 

=  / 

J  lb 

=  / 

J\\y 


aTQy 

l<t  brQyny 
aT[bx  Q\}y 
|<t  bT[bl  Qi]y 
,Thl 


d„y 


a‘blyi+a‘  Q\y 


jdny.  (35) 


I<1  bTbly i  +bJ  Q\y 
=o 

The  expression  (35)  can  be  further  simplified  as  follows: 
aTbl 


f  ax 
i|h|[<i  bTx 


d„x  = 


k 

vL 


IMI<i 


dny 

aTQiy 


(36) 


dn-iydyi. 


'iMi<i  y i 

The  first  integral  on  the  right  of  (37)  evaluates  to  V„,  the 
volume  of  a  unit  sphere  in  R'\  whereas  the  second  integral 
vanishes  (it  is  improper  but  the  Cauchy  principal  value  is  0). 
Therefore, 

aTb 


C  aTx  a1  b 

J\\x\\<\  bTx  "  bTb V" 


(37) 


/|h||<l  bTx 

where  Vn  is  given  by  (33).  This  completes  the  proof.  ■ 
Next,  the  relative  participation  of  the  k- th  state  in  the  i-th 
mode  is  evaluated  using  Definition  1  under  the  assumption 
above  on  the  distribution  of  the  initial  condition  x°.  Before 
proceeding,  we  recall  the  relationship  between  x°  and  z°: 


=  Vzu 


E  rjz°J- 

j= 1 


(38) 


A.  Participation  in  a  mode  associated  with  a  real  eigenvalue 
To  determine  the  participation  of  the  A-th  state  in  a  real 
mode  associated  with  a  real  eigenvalue  A.,-,  we  substitute  x?  = 
JJ o 


r;  ,c-; in  wy- 


%kj  —  E 


=  E 


ijA 

z° 

■ 


rJ7 o ' 
=1  rkzj 


=  E 


'  71  J  ,0 
lkrkZi 


+ 


VP 

lkrk 


E 

j= h  m 


l[r{E 


(39) 


Note  that  the  first  term  in  (39)  coincides  with  pki,  the  origi¬ 
nal  participation  factors  formula.  We  will  find  that,  in  general, 
the  second  term  in  (39)  does  not  vanish.  This  is  true  even 
in  case  the  components  x^x®,  •  •  •  ,x®  representing  the  initial 
conditions  of  the  state  are  assumed  to  be  independent.  This  is 
due  to  the  fact  that  the  second  term  involves  the  components 
of  z°  (i.e.,  Zj,Z2, •  •  •  ,Z«)  which  need  not  be  independent  even 
under  the  assumption  that  the  x?  are  independent,  due  to  the 
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transformation  z°  =  V  lx°.  This  was  overlooked  in  [1],  leading 
to  the  incorrect  conclusion  there  that  the  second  term  in  (39) 
vanishes. 

We  now  use  Lemma  1  to  simplify  the  expression  (39)  for 
the  participation  factor  for  the  k-th  state  in  the  ;-th  (real)  mode: 


S  (*°)  = 


l\x\  +  l‘2x 2  H - 1-  l‘nx® 


Denote  a  :=  ...  fh)1  and  b  :=  ■  ■  ■  JhY  •  The 

pected  value  of  g(x°)  is 

£{•?(/)}  =  [U0U  g(x°)f(x°)  dx° 


£ {•?(/)}  =  ^kv» 

aTb 

bTb 


*ki  =  lkrk  +  £  Z< 


i\T 


k'k 


j= i,  m 


ifvy 


2 


Tiki  = 


C I'k ) 

/>'(/ or 

(4)2 


D2- 


This  expression  is  easily  obtained  from  Definition  1  and 
Lemma  1  as  follows.  Recall  Definition  1,  the  general 
averaging-based  definition  for  state-in-mode  participation  fac¬ 
tors  for  the  case  of  a  real  eigenvalue: 


(40) 


n  (  1 

**»■  =  zi-4  +  e  iiriE  Hr 

j= i,  j&  i  j 

Substituting  z°  =  /'x°  into  (40)  yields 

=  «+  E 

7=1,  ^  X  J 

=  Vk  +  E  4^  {,?(*0)}  (4i) 

y=b  j¥=i 

where  g(x°)  takes  the  form 

/jX?  T-  /  7X2  3 - h  InX^ 


■  /'  r0  ■ 

**  := 

Zi 


(47) 


(48) 


Substituting  z"  =  /'x  in  (47)  yields 

f  Z[x°  1  f  liekx°  } 

*»  =  EW\=E\l^r}- 

Denote 

a  :=  l[ek  =  l[[ 0  ...  0  1  0...  0]r, 
b  :=  (/')J. 

Using  Lemma  1  and  the  normalization  kVn  =  1,  (48)  reduces 
to 


(42) 


K/d  — 


rb  l[ek(l‘)T  (l[)2 


bTb  v{vy  /'(/ iy1 


(49) 


which  is  exactly  (46). 


=  ‘t,a  *°  <43) 


B.  Participation  in  a  mode  associated  with  a  complex  conju¬ 
gate  pair  of  eigenvalues 

To  determine  the  participation  factor  for  a  state  in  a  complex 
mode ,  i.e.,  a  mode  associated  with  a  complex  conjugate  pair 
of  nonreal  eigenvalues  Xi,  X*,  we  use  the  second  case  of  (27): 


Using  Lemma  1,  which  applies  since  b  is  real,  and  the 
normalization  kVn  =  1  from  (31),  this  integral  reduces  to 


%ki  —  E  • 


=  E 


Vl+ik)4 

z°  +  z°* 

*•1  1  l 

.  Re{z?} 


(44)  Substituting  z®  =  Z'x°  in  (50)  yields 


Substituting  (44)  into  (41)  yields  a  key  result  of  this  paper,  a 
new  formula  for  the  participation  factor  for  state  Xk  in  a  real 
mode: 


Tiki  =  E 


Re®- 


IK 


Re{/‘x0} 
Equation  (51)  can  be  rewritten  as 


=  E 


R  e{/jK 

Re{/'}x° 


(50) 


(51) 


(45) 


Hr  e{4}(e*)rx° 

Kki  j  Re{/'}x° 


Remark  2:  Under  the  initial  condition  uncertainty  assump¬ 
tion  based  on  which  (45)  was  obtained,  the  participation  factor 
for  the  ;-th  mode  in  the  k-th  state  equals  the  participation  factor 
for  the  k-th  state  in  the  ;-th  mode  (i.e..  Tiki  =  Pki)  if  the  left 
eigenvectors  of  the  system  matrix  A  are  mutually  orthogonal, 
i.e., 

l]{l')T  =  0,  for  j,i=  1,2, .. . ,n,  i^j. 

This  is  a  very  restrictive  case  (which  applies,  for  instance, 
when  the  system  matrix,  being  real,  is  symmetric). 

Next,  we  derive  the  following  expression  equivalent  to  (45) 


=  R  e{4}£ 


(ek) 


(52) 


Re{/'}x° 

Next,  we  obtain  formulas  analogous  to  (45)  and  (46)  above, 
but  now  giving  participation  factors  measuring  participation  of 
states  in  a  complex  mode. 

First,  to  determine  a  formula  analogous  to  (45),  substitute 
Z-J  =  llx°  and  x°k  =  £"=1  rJkz°j  =  £"=1  r/Wx0  and  apply  and 
Lemma  1  to  obtain 

=  e-Rc{i m 


Re{/'}x° 

/;  [  Re  i :^x"' 


\  Re{/!> 


(46) 


=  Re{4-}£' 


E"=l  rJkljx° 
Re{/'}x° 


(53) 
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Let  a  =  | YJj=irJkl ^  and  /?  =  (Re{/'})r.  Since  /?  is  real, 
we  can  invoke  Lemma  1  and  the  normalization  kVn  =  1  to 
reduce  (53)  to 


^Ici 


(Vj=lr{v)(  Re{/‘})r 
Re^A'}  Re{Z,}(Re{Z,})r 


This  formula  can  be  rewritten  as 


(54) 


„  . p  tji i  (rkli)  (Re{/'})r  ,  £  Re{Z'}r/,L(Re{/-'})r 

1  ^  Re{/'}(Re{/f})r  Re{/«-}(Re{/-'})r 

(55) 


which  is  the  desired  form  analogous  to  (45).  We  observe  that  if 
this  formula  is  applied  to  a  simple  real  eigenvalue  A,-,  implying 
that  the  associated  eigenvector  can  also  be  taken  as  real,  the 
formula  indeed  reduces  to  the  formula  (45)  that  was  derived 
for  the  case  of  a  real  mode.  Thus,  formula  (55)  provides  a 
general  expression  for  state-in-mode  participation  factors  for 
systems  without  any  restriction  on  the  eigenvalues  besides  that 
they  are  distinct. 

Next,  we  obtain  another  equivalent  formula  for  the  Tiki  for 
a  complex  mode,  but  this  time  in  a  form  analogous  to  the 
expression  (46)  derived  above  for  the  case  of  real  eigenvalues. 
We  use  Lemma  1  to  simplify  the  expression  (52).  Taking  a  = 
ek  and  b  =  ( Re { /' })7  in  Lemma  1  and  using  the  normalization 
kVn  =  1,  (52)  reduces  to 


ttki 


R  e{/£} 

Re® 


(e*)r(Re{/'})r 

Re{/'}(Re{/'})7’ 

Re^} 

Re{/'}(Re{/'})7’' 


Finally, 


(Re{/j.})2 

'  Re{/'}(Re{/'})7’' 


(56) 


This  expression  is  an  alternate  form  of  the  newly  proposed 
formula  (55)  for  participation  factors  measuring  participation 
of  a  state  Xk  in  a  mode.  Both  expressions  apply  for  a  mode 
associated  with  a  complex  conjugate  pair  of  eigenvalues  A,- 
and  A*,  as  well  as  for  a  simple  real  eigenvalue  A 

Remark  3:  Although  care  was  taken  with  respect  to  se¬ 
lecting  units  in  deriving  the  formulas  for  state-in-mode  par¬ 
ticipation  factors,  the  final  formulas  are  not  themselves  in¬ 
dependent  of  units.  The  independence  with  respect  to  state 
variable  units  that  occurs  in  the  definitions  of  mode-in-state 
participation  factors  is  a  fortunate  coincidence  for  certain 
choices  of  initial  conditions  or  under  certain  initial  condition 
symmetry  assumptions.  However,  no  such  coincidence  occurs 
in  the  quantification  of  state-in-mode  participation.  One  way 
of  viewing  this  is  as  follows.  Units  are  important,  in  the  sense 
that  they  should  be  chosen  so  that  a  unit  variation  in  the  initial 
condition  of  any  state  variable  has  a  similar  likelihood  as  a  unit 
variation  in  the  initial  condition  of  any  other  state  variable  of 
the  system.  That  is  the  spirit  of  the  assumption  made  above 
that  the  distribution  of  initial  conditions  of  the  state  vector  can 


be  mapped  by  changes  of  units  to  a  uniform  distribution  on  a 
unit  sphere  in  R". 

Next,  we  revisit  Examples  1  and  2  using  the  newly  derived 
formula  for  state  in  mode  participation  factors,  and  compare 
the  results  to  the  participation  factors  obtained  using  the 
original  definitions.  Note  that  all  eigenvalues  in  these  examples 
are  real,  so  the  formula  (45)  applies  as  a  (new)  measure  of 
state-in-mode  participation  factors  Tiki  (as  does  the  equivalent 
formula  (46)). 

Example  1  Revisited 

For  Example  1,  the  participation  factors  for  states  x\  and  X2 
in  mode  1  based  on  the  new  formula  (45)  are 

(i a  —  d )2  b 2 

KU  =  (a-df  +  b2  “  %21=  (a-d)2  +  b2’ 

respectively.  The  participation  factors  for  states  x\  and  X2  in 
mode  1  based  on  the  original  formula  are  p\\  =  1  and  pi \  =  0, 
respectively. 

As  we  observed  previously  in  our  discussion  of  Example 
1,  the  original  formula  for  participation  factors  erroneously 
indicates  that  the  participation  of  state  X2  in  mode  1  is 
zero.  The  coupling  between  state  X2  and  state  x\  in  the 
system  dynamics  is  not  reflected  in  the  original  formula  for 
participation  factors  (the  pki),  whereas  this  coupling  between 
state  variables  is  reflected  in  the  result  of  applying  the  new 
formula  (for  the  7t *,■). 


Example  2  Revisited 

For  Example  2,  the  participation  factors  for  states  x\  and  X2 
in  mode  2  based  on  the  new  formula  (45)  are 

1  ,  1 

7tl2=2  and  n22  =  ^, 

respectively.  The  participation  factors  for  states  x\  and  X2  in 
mode  2  based  on  the  original  formula  are 


P 12  = 


1 

1  -d 


and 


- d 

P22  =  T^d' 


respectively. 

The  results  using  the  new  formula  more  faithfully  reflect  the 
relative  contributions  of  the  initial  conditions  of  the  two  state 
variables  to  the  evolution  of  mode  2,  which  is  given  explicitly 
by  the  formula 

22 0 )  =  yTj  (x'l+xtye12*,  (57) 

Here  it  is  clear  that  z.2(t)  is  equally  influenced  by  x®  and  x® 
since  it  depends  on  the  initial  condition  x°  through  the  sum 

Xj  +Xi. 


VI.  A  Numerical  Example:  A  Two-Mass 
Mechanical  System 

Consider  the  translational  mechanical  system  depicted  in  Fig¬ 
ure  1,  where  y\(t)  and  y2{t)  denote  the  displacements  of  mass 
1  and  mass  2,  respectively,  from  the  static  equilibrium  [13]. 
The  system  parameters  are  the  masses  m\  and  m2,  viscous 


9 


damping  coefficients  c i  and  C2 ,  and  the  spring  constants  k\ 
and  k 2. 

A  state  space  representation  is  obtained  by  defining  the 
system  states  as 

*1  (0  =  yi{t) 
xi  (t)  =  y  i=i'i 
+s(0  =  y2(t) 
x4(f)  =  y2=xi- 


The  system  dynamics  is  described  by  the  linear  time-invariant 
differential  equation  [13] 


x  =  Ax 

where 

1  0  0' 

C\  +C2  j£2_  £l_ 

m\  m\  m\ 

0  0  1' 

c 2  h  £2 

mi  »>2  m2  - 

With  the  system  parameters  selected  as  m\  =  39  kg,  mi  =  17 
kg,  ci  =  19  Ns/m,  C2  =  33  Ns/m,  k\  =  374  N/m,  and  &2  =  196 
N/m,  the  system  state  dynamics  matrix  A  becomes 


A  = 


0 

k  I  ki 
m  ] 

0 

±L 

mi 


A  = 


0  1 

-14.6154  -1.3333 

0  0 

11.5294  1.9412 


0  0 

5.0256  0.8462 

0  1 

-11.5294  -1.9412 


respectively.  Note  that  since  X2  =  X\  and  X4  =  X3,  we  have 
that  r2  =  rl*  and  r4  =  r3*,  where  an  asterisk  denotes  complex 
conjugation.  The  left  (row)  eigenvectors  of  the  system  matrix 
A  associated  with  Xi  and  X3  are 


l1 


/3 


-0.3122  +  ;1.2059 
0.4806  +  y'0.0539 
0.2760  +  y'0.7884 
0.3730 -;'0.0171 

0.3134 -;2.4251 
-0.4824 -./'0. 1776 
-0.2771 +;1. 3357 
0.2530  +  y'0.1903 


respectively,  and  l2  =  l1*;  /4  =  /3*. 

The  magnitudes  of  the  original  participation  factors  pia 
evaluated  using  (23)  for  this  example  are  given  in  Table  I. 
The  state  in  mode  participation  factors  evaluated  using  the 
new  formula  (56)  are  given  in  Table  II. 


TABLE  I 

Participation  factors,  pki,  based  on  original  formula. 


mode  1 

mode  2 

*1 

\pn  |=0.2420 

|pi2|  =0.3050 

*2 

\P2i  j  —  0.2184 

|p22|  =0.2900 

*3 

|P31 1  =0.2874 

\P32 1  =  0.2404 

X4 

|P4i  |=0.2987 

P42I  =  0.2523 

TABLE  II 

State  in  mode  participation  factors,  nki,  based  on  new  formula. 


mode  1 

mode  2 

X\ 

Jin  =0.1792 

7ii2  —  0.2082 

*2 

7I2i  =  0.4248 

7122  =  0.4934 

x3 

JI31  =  0.1401 

jl32  =  0.1628 

X4 

H41  =  0.2558 

7142  =  0.1357 

Fig.  1.  Mechanical  system. 


The  eigenvalues  of  A  are  X\i  =  — 0.217  ±  /2. 315  and  A-3,4  = 
—  1.4203  +  j4.2935.  We  denote  the  modes  associated  with 
the  eigenvalues  as  21(f)  and  22(1)-  respectively.  The  dominant 
mode  is  the  one  associated  with  the  complex  conjugate  pair 
of  eigenvalues  closest  to  the  imaginary  axis,  i.e.,  X\p..  Modes 
associated  with  eigenvalues  close  to  the  imaginary  axis  are  of 
considerable  interest  as  they  can  be  used  as  an  indication  of 
closeness  to  system  instability.  Therefore,  our  emphasis  in  this 
example  will  be  on  21(f). 

The  right  (column)  eigenvectors  of  the  system  matrix  A 
associated  with  A-i  and  X3  are 


r 


1 


r 


3 


0.0124  -JO.  1938 
0.4461  +  J0.0708 
-0.0321  —./0.3425 
0.7999 

0.0716  +  7O.IO2I 
-0.5402  +  JO.  1624 
-0.0553  -JO.  1673  ’ 

0.7970 


(58) 


We  observe  that  the  participation  factors  given  in  Table  I, 
which  are  calculated  using  the  original  definition  of  partic¬ 
ipation  factors,  differ  from  the  state  in  mode  participation 
factors  in  Table  II  calculated  using  the  new  formula.  For 
instance,  according  to  Table  I,  the  state  that  participates  most 
in  mode  1  is  X4,  whereas  according  to  Table  II,  the  state 
that  participates  most  in  mode  1  is  X2-  To  demonstrate  that 
state  X2  participates  more  than  other  state  variables  in  mode 
1,  we  calculate  the  evolution  of  mode  1  (21(f))  due  to  different 
settings  in  the  initial  conditions.  Specifically,  21  (f )  is  calculated 
for  the  following  set  of  initial  conditions:  x°  =  [0.1,0,0,0]r, 
x°  =  [0,0.1,0,0]r,  x°  =  [0,0,0.1,0]r,  and  x°  =  [0,0,0,0.1]r. 
We  select  the  initial  conditions  in  this  way  (0  in  all  but  one  of 
the  state  variables)  to  distinguish  the  influence  of  each  of  the 
state  variables  on  mode  1 .  The  simulation  results  are  depicted 
in  Figure  2,  which  shows  plots  of  Re{2i(f)}  for  the  various 
initial  conditions.  Figure  2  shows  that  the  initial  condition 
component  x®  gives  the  largest  effect  on  mode  1  at  f  =  0 
compared  to  all  other  state  variables.  This  agrees  with  what  is 
predicted  using  the  new  formula  for  state  in  mode  participation 
factors  (see  Table  II).  In  other  words,  mode  1  can  be  excited 
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by  an  initial  condition  on  any  of  the  state  variables,  however, 
the  highest  excitation  at  t  =  0  comes  from  x®  and  the  least 
excitation  comes  from  *3. 

The  state  in  mode  participation  factors  calculated  based  on 
the  new  formula  derived  in  this  paper  can  be  used  to  determine 
the  relative  degrees  by  which  system  states  excite  a  particular 
mode  in  the  system.  This  can  be  useful  in  stability  monitoring 
applications.  For  example,  in  stressed  electric  power  systems, 
typically  there  is  one  critical  mode  that  needs  to  be  monitored. 
If  the  system  is  operating  exactly  at  an  equilibrium  point  and 
not  influenced  by  disturbances,  this  mode  cannot  be  observed 
in  the  outputs  of  the  system.  In  order  to  monitor  the  critical 
mode  and  detect  closeness  to  instability,  a  small  perturbation 
signal  is  applied  to  the  system  and  the  response  is  measured. 
The  state  in  mode  participation  factors  can  help  in  selecting 
a  location  for  applying  the  perturbation  signal  in  order  to 
achieve  the  highest  excitation  in  the  critical  mode  to  achieve 
the  clearest  possible  indication  on  how  close  the  system  is  to 
instability. 


Consider  the  general  expression  (22)  for  the  mode  in  state 
participation  factor  p repeated  here  for  convenience: 

n  f  x0  j 

Pki  =  l'A+  £  iykE  .  (61) 

m  [xk ) 


A  typical  term  in  the  summation  on  the  right  side  of  (61)  takes 
the  form 


Ml 

E  \  of 

I  k>j±k 


[  4/(x°)  dx° 

vn<i*r 

k[ 

Jm\<i 


lh°ll<i  4 


Denote 


a  :  =  (0,...,0,^l^,0,...,0y, 

jth 

b  :=  (0,...  ,0,^,0,.  ..,0)r. 

k‘h 


(62) 


The  integral  in  (62)  can  be  expressed  as 


v-0 


=  k 


L 


T  0 
a  x 


|  |.y°  1 1  <  1  bTx° 


Using  Lemma  1,  this  expression  reduces  to 

aTb 


bTb 


dx 


=  0. 


(63) 


(64) 


Therefore,  the  second  term  in  (61)  vanishes  and,  under  the 
new  assumptions,  the  mode  in  state  participation  factors  are 
still  given  by 


Fig.  2.  The  effect  on  mode  1  of  a  perturbation  of  0. 1  away  from  equilibrium 
in  the  initial  condition  component  x ®  (solid),  x\  (dotted),  a)  (dash-dot)  and 
*4  (dashed). 

VII.  A  Further  Remark  on  Mode  in  State 
Participation  Factors 

In  this  section,  an  additional  result  is  given  that  relates 
only  to  mode  in  state  participation  factors.  The  result  is  that, 
in  the  averaging  formulation,  an  additional  possible  set  of 
initial  condition  uncertainty  assumption  is  found  to  also  lead  to 
the  traditional  participation  factors  formula  for  mode  in  state 
participation  factors. 

In  Section  IV,  we  showed  that  when  the  components  of  the 
initial  condition  vector  x^jx®, . . .  ,x®,  are  independent  random 
variables  with  zero  mean,  the  participation  of  mode  i  in  state 
Xk  is  given  by 

Pki  =  1‘A-  (60) 

In  this  section,  we  show  that  this  expression  remains  valid  if 
the  components  of  the  initial  condition  x°,  x(j.  j  =  1,2, .. .  ,n, 
are  assumed  to  be  jointly  uniformly  distributed  over  the  unit 
sphere  in  Rn  (see  (28)  for  the  expression  of  the  probability 
density  function  in  this  case).  Note  that  under  this  assumption, 
the  random  variables  xj  are  no  longer  independent. 


Pki  =  A  (65) 

We  can  therefore  conclude  that  (65)  is  a  valid  formula  for 
mode  in  state  participation  factors  under  any  of  the  following 
assumptions  on  the  initial  conditions: 

1)  The  initial  condition  x°  is  taken  to  lie  in  an  uncertainty 
set  §  which  is  symmetric  with  respect  to  each  of  the 
hyperplanes  x?  =  0,  k  =  1,2, ...  ,n. 

2)  The  initial  condition  components  are  independent  ran¬ 
dom  variables  with  marginal  density  functions  which  are 
symmetric  with  respect  to  x?  =  0,  k  =  1,2, . . .  ,n 

3)  The  initial  condition  components,  x(],  j  =  1,2 are 
jointly  uniformly  distributed  over  a  sphere  centered  at 
the  origin. 

VIII.  Conclusions 

We  have  presented  a  new  fundamental  approach  to  quan¬ 
tifying  the  participation  of  states  in  modes  for  linear  time- 
invariant  systems.  We  have  proposed  that  a  new  definition 
and  formula  replace  the  commonly  used  participation  factors 
formula  for  measuring  participation  of  states  in  modes,  while 
recommending  the  previously  used  participation  factors  for¬ 
mula  be  retained  as  a  quantification  of  participation  of  modes 
in  states.  The  analysis  presented  in  this  paper  uses  averaging 
over  an  uncertain  set  of  system  initial  conditions.  The  analysis 


11 


led  to  the  conclusion  that  participation  of  modes  in  states 
and  participation  of  states  in  modes  should  not  be  viewed 
as  equivalent  or  interchangeable.  Examples  were  given  to 
demonstrate  that  the  original  formula  for  participation  factors 
is  not  convincing  as  a  measure  of  state  in  mode  participation. 
Moreover,  these  examples  demonstrated  the  applicability  and 
usefulness  of  the  new  formula  for  state  in  mode  participation 
factors. 

It  is  interesting  that  while  the  problem  addressed  in  this 
paper  relates  to  a  very  simple  and  well  studied  class  of 
systems,  considerable  effort  was  required  to  revisit  what  may 
seem  to  be  a  basic  issue,  namely  modal  participation.  Indeed, 
it  appears  that  work  is  needed  on  further  related  matters,  such 
as  implications  of  the  results  for  sensor  and  actuator  place¬ 
ment,  for  system  monitoring  to  detect  impending  instability, 
for  order  reduction,  and  possibly  for  coherency  studies  of 
power  networks,  among  other  issues.  Also,  from  the  results 
in  Figure  2  for  the  mechanical  example,  the  relative  sizes 
of  modal  participations  at  the  initial  time  instant  might  differ 
from  those  over  a  time  interval.  For  some  applications,  it  may 
be  desirable  to  have  analytical  measures  of  modal  participation 
that  quantify  modal  participation  over  time. 
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